THE GOURSAT PROBLEM FOR THE EINSTEIN- VLASOV 
SYSTEM: (I) THE INITIAL DATA CONSTRAINTS 



CALVIN TADMON 

^ ■ DEPARTMENT OF MATHEMATICS AND COMPUTER SCIENCE, FACULTY OF 

^ ! SCIENCE, UNIVERSITY OF DSCHANG, P. O. BOX 67, DSCHANG, CAMEROON 

(N : 

Q^' Abstract. We show how to assign, on two intersecting nuh hypersur- 

D . faces, initial data for the Einstein- Vlasov system in harmonic coordinates. 

' As all the components of the metric appear in each component of the 

. stress-energy tensor, the hierarchical method of Rendall can not apply 

CO ' strictly speaking. To overcome this difhculty, an additional assumption 

. have been imposed to the metric on the initial hypersurfaces. Conse- 

^ ' quently, the distribution function is constrained to satisfy some integral 

O^' equations on the initial hypersurfaces. 

^ ! 1. Introduction 

^ I This work is devoted to the resolution of the constraints problem associ- 

OO ■ ated to the characteristic Einstein- Vlasov (EV) system on two intersecting 

^ . null hypersurfaces. The interests and physical motivations for studying such 

g ; problems have been widely mentioned in[21lll[SlEl[ini[I21[ISl[Il[IZllIB]. It 

^ I is well known that the EV system is not an evolution system as it stands. In 

^ ' order to obtain a hyperbolic system, one needs to impose some supplemen- 

>• . tary conditions called gauge conditions which, due to the deep structure of 

^ i the system, must satisfy the following properties: 

I {{) whenever these gauge conditions are fulfilled everywhere in the space- 

time, the EV system reduces to a non-linear hyperbolic system called the 
evolution system. 

(ii) whenever the associated evolution system is satisfied everywhere in the 
space-time and the gauge conditions are satisfied on the null hypersurfaces 
that carry the initial data, then these gauge conditions and the complete EV 
system are satisfied everywhere. 

It therefore follows that when the choice of gauge conditions is made, the 
initial value problem for the EV system is naturally decomposed into two 
parts called the evolution problem and the constraints problem. 
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The resolution of the evolution problem is equivalent to the resolution 
of the reduced non-linear hyperbolic system obtained from the EV system 
thanks to the choice of the gauge conditions. Due to the gauge conditions 
the data for the reduced EV system can not be given freely. It is necessary 
to construct, from arbitrary choice of some components of the gravitational 
potentials (called free data) on the initial null hypersurfaces, all the initial 
data such that the solution of the reduced EV system with those initial data 
satisfies the gauge conditions on the initial null hypersurfaces. The construc- 
tion of such data is referred to as the resolution of the constraints problem. 
Through all the work we will use harmonic gauge for the gravitational field. 

We now proceed to survey some relatively recent works known about char- 
acteristic initial value problems with initial data prescribed on two intersect- 
ing null hypersurfaces often referred to as the Goursat problems. In 1990, 
A. D. Rendall p/Tj published a C°° existence and uniqueness result for quasi- 
linear hyperbolic systems of second order with data prescribed on two 
intersecting null hypersurfaces. Using the harmonic gauge, the author ap- 
plied the result obtained in [T7] to solve the characteristic initial value 
problem for the Einstein equations in vacuum and with relativistic perfect 
fluid source. For sake of more physical applications, it is known that, for 
Partial Differential Equations (PDE), solutions of finite differentiability or- 
der are more important than those of infinite differentiability order. In [T7j 
section 7 the author mentioned briefly how results of finite differentiability 
order can be obtained for data of finite differentiability order although proofs 
were not given. In 1990, H. Miiller zum Hagen [H] used Sobolev type in- 
equalities to derive energy inequalities that enable him solve, in weighted 
Sobolev space (results of finite differentiability order), the characteristic ini- 
tial value problem for linear hyperbolic systems of second order. He also 
predicted an existence and uniqueness result for the quasilinear case. Apart 
from the fundamental papers [H] and [17], some other works on character- 
istic initial value problems with initial data prescribed on two intersecting 
null hypersurfaces can be found in [21 [3|, El 13 [101 [121 [H]. As pointed out 
by H. Andreasson [1], A. D. Rendall [17] and M. Fjallborg [11], unlike some 
known models, the Einstein- Vlasov model has a very nice feature in General 
Relativity and Kinetic Theory since the stress-energy tensor fulfills, with- 
out any supplementary assumption, all the physical necessary energy condi- 
tions i.e. the weak energy condition, the dominant energy condition and the 
strong energy condition as well as the non-negative sum pressures condition. 
This situation, coupled with the importance of characteristic initial value 
problems mentioned at the beginning, motivates us to study the constraints 
problem associated to the characteristic EV system. When attempting to 
solve the constraints problem for the characteristic EV system by the hier- 
archical method of Rendall (see [T71 [H]), a crucial obstacle occurs due to 
the complicated form of each component of the stress-energy tensor where 
all the components of the metric to be constructed appear. The novelty of 
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our work resides in the fact that we have worked out this difficulty through a 
supplementary judicious assumption imposed to the gravitational potentials 
on the initial hypersurfaces. As a consequence of the additional assumption 
on the gravitational potentials, the distribution function can not be given as 
free data, it must satisfy some integral equations. Another advantage of this 
paper is that, unlike the work of Rendall [HI [H], many delicate calculations 
and expressions are given in details in such a way that we can foresee promis- 
ing resolution of the global characteristic EV system using for example tools 
that are similar to those of G. Caciotta and F. Nicolo [31 S]. To reduce the 
length of the paper, the evolution problem for the characteristic EV system 
is out of the scope of the present work and will be solved in a forthcoming 
paper. The paper is organized as follows. In section 2, we give some pre- 
liminaries about the EV system. The complete form as well as the reduced 
form (in harmonic coordinates) of the EV system are written. A new form 
of EV system is derived with appropriate unknowns an variables. This new 
form of EV system is suitable for the resolution of the constraints problem. 
The concern of section 3 is the resolution of the constraints problem for the 
characteristic EV system i.e. the construction of the initial data for the re- 
duced EV system such that the harmonic gauge conditions are satisfied on 
the initial null hypersurfaces. For sake of simplicity and clarity, only the 
case of C°° data will be discussed. Data of finite differentiability order may 
be constructed in Sobolev type spaces using energy inequalities and other 
classical tools as described in [9l [TOl [HI [17] and references therein. An ap- 
pendix D is provided at the end of the work and is devoted to the treatment 
of the constraints integral equations which must be satisfied by the distribu- 
tion function. It would be of interest to investigate whether the additional 
assumption on free data as well as the constraints integral equations have a 
particular physical meaning. 



2. The Einstein- Vlasov (EV) system 

2.1. The complete form of the EV system. The geometric framework 
is a four dimensional differentiable manifold Ai, endowed with a hyperbolic 
metric ^of signature — +++. The manifold (A^, ^) is called a space-time. Ai 
is assumed to be orientable and of class C°°. Throughout the remainder of 
the work, commas will be used to denote partial derivatives e.g. ^ij^k = ^f- 
Roman indices i,j, ... run from 1 to 4 while Greek ones a, /3, ... run from 3 to 
4. Einstein convention on repeated indices is used i.e AiB^ = '^AiB\ The 

i 

Einstein- Vlasov system is written as follows (see [H [71 [Til US]) 
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where are the covariant components of the metric ^. They constitute the 
unknowns for the Einstein equations. Rij are the covariant components of 
the Ricci tensor and R is the scalar curvature of the metric ^. In the local 
coordinates (y*) they are given as follows 

Rij = Rikj = ^ij,k ~ ^ik,j + ^kl^ij ~ ^^i^lk^ — 9^ Rij^ (2-2) 
where T^j are the Christoffel symbols of the metric ^ i.e., 

r}fc = {gik,j + gij,k - gjk,i) , (2.3) 

are the contravariant components of ^ i.e., 



/ is the distribution function (or the particle number density function) which 
constitutes the unknown for the Vlasov equation. / is a non-negative real 
valued function defined on F {Ai ) , where 

F [M) = {q = [q') G T, : g,, (y) q^ = -m\ < g^} , (2.5) 

with Ty = TyAi. The Vlasov equation symbolizes the conservation of the 
number of particles along the trajectories across the hypersurfaces of F {M.) 
in the case where there is no collision between particles. Tjj are the covariant 
components of the stress-energy (or energy-momentum) tensor which is the 
source of the gravitational field created by the particles. In contravariant 
components the stress-energy is defined by the following relation (see [7]) 

^ 1 

r'{y) = - [ f{y,q)qW^-^d\ (2.6) 

JFy Ql 

where Fy = {q = (g*) G Ty : gij (y) q'q^ = -m^, < g^}, d^q = dq^ A dq^ A 
dq'^, 1^1 is the modulus of the determinant of (jj, 



2.2. The reduced EV system. The Einstein equations as they stand are 
not hyperbolic but in harmonic coordinates they read (see |7j) 



where 



Ri = (2.7) 

R-ij = Rij ~ 2 {9i''^d + ^jk^'^i) = -i^'^gij^mk + Qij- (2.8) 

Here Qij is a rational function depending on the metric components and 
their first order derivatives (see [ini 120]). 

= g^^n. (2.9) 
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So the reduced EV system reads 

9ij,mk + Qij = Tij, , , 

where 

= -%q^q'. (2.11) 

2.3. Appropriate unknowns and variables. As a relativistic speed is 
bounded, we think that it is convenient to choose on the mass shell, coordi- 
nates with bounded domain (see also ff\). Let y e U, U is the domain of a 

A 

local chart in Ai. Set w = A = 2,3,4 and denote by My the image in 
of Fy by the mapping (g*) t— )■ (w^). Assume the following hyperbolicity 
conditions on (jjij) . 

Assumption h: The metric {'Qij) is uniformly hyperbolic and the hyper- 
surfaces = Cte are uniformly spatial i.e. 

3a,be{0,oo):a^\^f<gAB^^e<b'\^f where l^l' = E {^^)\ 

A=2 

—Qii > «^ and — g^^ > c? . 

(2.12) 

Proposition 1. ii) Under assumption ( \2.1S^ . the Vlasov equation reads 

?'I4 + Q^I4 = 0. (2.13) 



(ii) Under assumption {\2.1<^ . My is a bounded domain in M such that 
My C M, where M is a fixed compact domain in M^. The stress-energy 
tensor is given as follows 

f''{y) = —J fiy,w)q'q^{q')'\gr^d'w, (2.14) 

J My 

where d^w = dw"^ A dw'^ A dw'^ , f {y,w) is the expression of f {y,q) in the 
local coordinates {y,w). 

Proof. See [7]. □ 



Remark 1. In the expression ^2.1^^ of the stress-energy tensor, we would 
like to write f {y,w) q^q^ {o^)'^ as if {y,w)w'^\ where w^^ does not depend on 
(gij). To do so, we proceed to the following change of the unknown distri- 
bution function by setting f(y,w) = (p {y,w) {q^) ^. So we must have w^^ 

~ W' 

Proposition 2. Under the change f {y,w) = (p{y,w){q^) ^, the stress- 
energy tensor is given as follows 
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r'{y) = ^l ^{y,w)w'^\g\U'w, (2.15) 

'My 



where 

w 

The Vlasov equation becomes 



Proof. See m. □ 

Remark 2. The expression i \2.15\i of the stress- energy tensor is not appro- 
priate since the domain My depends on y and makes it difficult to differen- 
tiate T^^ even in the distributional sense. It appears therefore judicious to 
transform this domain in order to make it independent of y. 

Assume the following decomposition of 'g. 

Assumption h': The spatial part of {c]ab) is decomposed as follows 



SabY^Y^ = Y,{\'iY^)\ (2.17) 

B=2 

where are functions that depends smoothly (C°° for instance) on the 
components qab of the metric. Set 

v^={-rf'\'A[^^ + 9'^]. (2.18) 

with 

g'^ = -^g'^. (2.19) 

fl2.18p is equivalent to the following relation 

w^={Xi)-\,^-rr'^-g'^ (2.20) 

where (A;^) ^ are the components of the inverse of the matrix (A j) . 

Proposition 3. The image of My by the mapping (w"^) i— )■ [v^^ is the unit 
open ball B in . In the parameters (f"^), the energy-momentum tensor 
reads 

T'Hy) = \ I ^{yMv''\W {-f'T'" \g\~'^d\, (2.21) 

Jb 

where d^v = dv'^ Adv^ Adv"^, v^^ = ■^^j Ifi'l is the modulus of the determinant 
of ((Jab), fiy^v) is the expression of(f{y,w) in the local coordinates {y,v). 
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The Vlasov equation reads as follows 



di£_ I 
9yi + 



lA 



Proof. See f^. 



w 



dip 



6(gi)"'gv = o. 



(2.22) 



□ 



Remark 3. Y. Choquet-Bruhat [7] used assumption [2.11] and a variant 
of assumption i \2.12i) to treat the ordinary Cauchy problem for the EV sys- 
tem. But in the characteristic case, those assumptions are not appropriate 
and they need to be recast. We proceed to the desired adaptation through a 
judicious change of local events variables (y*). 

Proposition 4. Let (y*) be a local coordinates system on Ai in which the 
components (jjij) of the metric satisfy assumption i\2.12^ . Set 



X 

x" 



1 (..I 



2 ' y ) . 

y^, a = 3,4. 



X 



iv' - y') 



1 (n.l 
2 



In the coordinates system {x,p), the EV system reads 

idf^ I pidf ^ g 



(2.23) 



(2.24) 



where 



dx^ dxi 3kl {y) 5 R-ij (^) gti QxJ Rkl iv) i 

^1 fin 



9ij [X) Q^i Q^j : 
R (x) = g'^ (x) Rij (x) , Tij (x) 



ai-i '^ki ^y) 



P — Qykq J -r — Qyk^ 

r!j, = 2^'*'^ {9ck,j + Qcj.k - 



-T)j^p^p^, 



jk — 2^ Vycfej -r ycj,k — 9jk,c) ■ 

Proof. A direct calculation leads to the desired equations. 



(2.25) 



□ 



Remark 4. The change of local coordinates i\2.23\) preserves the harmonic- 
ity. In other words, we have the following equivalence 



Vz = 1, 2, 3, 4, r = = o U» (Vz = 1, 2, 3, 4, V = g^'Vl^ 



In fact, the following relations hold 

fi = ri + r2, 
f " = r", a 



r2 = ri 

3,4. 



0). 

(2.26) 

(2.27) 



Proposition 5. The stress-energy tensor {\2.21\] is given in the local coordi- 
nates (x*) and the local parameters {v^) as follows 



rpl 



\x 



2m? 



(2.28) 
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where v^^ = = ! \9\ ^he modulus of the determinant of (gij), 

\g\ is the modulus of the determinant of {c]ab), f{x,v) is the expression of 
ip{y,v) in the local coordinates {x,v). The Vlasov equation (\2.23t} becomes 

with 

= (Xi)-' i-g'T' - 9'^, (y) = 9'' (x) + 2g'^ (x) + g'' (x) , 
H^{x,v) = + H^x,v) = l{l-w^), H''{x,v)=w", a = 3,4, 

LC (a;, v) = i i-g^'r & v) - \ (-?")^ A^/" (x, v) , F (x, v) = \w^F. 

(2.30) 

Here 

F = (r?, + r}0 (i + ^^)V(rL + ry (i-^^)V2(rL + ry (i-M^ 
+4 (rL + r\,) (1 + w^)+A {rl, + rl,) (i -w^)+a {rl + rij 



and 
with 



f {x,v) =ll + wHl r{x,v) =1^ -w'^ll a = 3, 4, 
+2 (rL - T\,) (1 - {w^Y) + 4 {Tl, - T\,) (1 + w^) 



+2 (rL + r\,) (1 - {w^Y) + 4 {rl, + r},) (i + w') 
+4 (r^, + r^,) (1 - + 4 (r^^ + rij 

+2r?2 (1 - M') + 4r?A^Mi + 

+4r^;,w^ (1 - u;^) + AT'^^w^w''. 

Proof. It is straightforward though lengthy. □ 

Remark 5. The EV system i \2.24^ in the local coordinates (x, f) reads as 
follows 

R,,-^Rg,,=T,,, W^ + L^^ + F^ = 0. (2.31) 
The reduced EV system i\2.10^ in the local coordinates {x,v) reads as follows 
Rij = T,„ H^^ + L^'^ + F^ = 0, (2.32) 



INITIAL DATA FOR THE CHARACTERISTIC EINSTEIN- VLASOV SYSTEM 9 



where ^ ^ 

Rij = Rij — - {gki^j + dkj^'^i} = '~-^3^"^9ij,mk + Qij- (2.33) 

3. The constraints problem for the EV system 

The task here is the construction of initial data for the reduced EV such 
that the constraints F'^ = are satisfied on G^UG'^, where G^ and G^ are the 
hypersurfaces in defined by = and = respectively. We will need 
G^ = {x e G'^ : < + < T} , T > 0, u = 1,2. Here the problem is 
much more difficult than in [Tni[I7|. This difficulty has something to do with 
the appearance of all the components of the metric in any component of the 
stress-energy tensor. To overcome this toughness, we add a supplementary 
assumption on the metric along the initial hypersurfaces. All the same, we 
try to mimic, as far as possible, the hierarchical method of Rendall fTfi 
[18] . We will see that the supplementary assumption has as consequence to 
force the distribution function to satisfy specific integral equations on the 
initial hypersurfaces. The resolution of the integral equations derived can 
be achieved under suitable conditions (see appendix D). Nevertheless it is 
still to be investigated whether our additional assumption has a particular 
physical meaning. The construction will be made in a standard harmonic 
coordinates system. The existence of such standard harmonic coordinates 
system has been established by A. D. Rendall |17j . 

Let us now adapt the method of Rendall to construct C°° initial data 
for the EV system. The assumptions under which the work is achieved are 
such that only the data gap and gi2 have to be constructed on U G^, the 
relations F'^ = have to be arranged on G^ U G^, the relations (722,1 = 2(712,2 
and gii^2 = 2(7i2,i have to be established on G^ and G^ respectively. The 
construction of the data is done fully on G^ and it will be clear that data on 
G^ are constructed in quite a similar way. The first level of the hierarchy is 
now described. 



a matrix with determinant 1 



3.1. Construction of gap and g^ on G^, relations F^ = and (722,1 

^33 ^34 

at each point. Set gap = ^hap, where f2 > is an unknown function called 
the conformity factor. Assume as in [TOj [17] that 



2fi'i2,2- Let T e (0, 00), {ha^ 



g22 = g23 = g24 = o on G^. (3.1) 

The additional assumption on free data is the following 

5'ii = fi'13 = 5-14 = on Gy. (3.2) 



On G^ it holds 



gi2g'^ = l, g''=g'" = 0, 

g'' = g"' = 0, gxpg^' = Sl ^^"^^ 
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3.1.1. Expression of R22 andT22- 
Proposition 6. On G\, it holds that 

R22 = W^g''^gap,2 (2^12,2 - 522,1) + \gfgxp,2 - \ (5"^5a/3,2) 



^22 = (512) K22, 



(3.4) 



where 



i^22(a;) = ^ / v{x,v){l + v^fd^v. (3.5) 



Proof. See appendix A. □ 

Assume 521 = for i 7^ 1 on G\, 522,1 = 2512,2 on G\,. Then = is 
equivalent to (see [T0| [T7]) 



512,2 = 2^12-^- ^^-^l 



The equation 



:gfgo.p,2 - - (5"^5«A2) 2 = T22, (3.7) 



la/3 1 
^5,2 fi'a/3,2 - 2 

provides the following non linear second order ODE with the conformity 
factor VL as unknown 

^) ' + \K,,h:i - 2 (^) ^ = 2i^22 (512)^ . (3.8) 

If we set $7 = e^, then the following system of ODE is derived from (13. 6p 
and (13. 8p in order to determine the conformity factor together with 512 

2^22 = - {¥2? - 2K22 (512)' + \Kp,2hf, ^3 
5l2,2 = 29l2y,2- 

Let T G (0, 00). Assume h^,s, K22 E C°° (G^) . Take Vq, Vi, Wq E (E) , 
where E = G]. fl G^. Then there exists Ti G ]0,T] such that (13. 9p has a 
unique solution (^,512) G (G^^J x (G^J satisfying V = Vq, V2 = Vi, 
5i2 = on E. This follows from known local existence and uniqueness 
results concerning non-linear ODE with G°° data in Banach spaces. 

3.1.2. The condition 522,1 — 2512,2 = on G^^. On G^^, the reduced equation 
-R22 = ^22 is equivalent to the following homogenous ODE with unknown 
522,1-2512,2 (see 

(5^^)^512,2 (522,1 - 2512,2) - 5^^ (522,1 - 2512,2)^2 = 0- (3-10) 

Assume 522,1 — 2512,2 = on E. Then 522,1 — 2512,2 = on G^^ and so E^ = 
on Gy^. 
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3.2. Relations = on G^^. We seek for a combination between 
and that will provide an homogenous ODE on with unknown F". 

Proposition 7. On G^, it holds that 

R2a + ^gafiT^2 + (^9^^ 9 12,2 9 a 13 + ^9al3,2^ = V'a, (3.11) 

(3.12) 

where 

i^a = 1 i9'^) 912,2 [-2g'^gi2,a + 9''' (2(?a^,e - 9t.B,a)\ 

+ \g^p,2 [-2gP^g^^gi2,x + g^^g^' {2gx^,e - 9,.e,x)] 

+ \ {9^^9a^,2) , - 3 {9''gi2,2)^^ - 1 {9''? 9l2,29l2,a 

y^"^ {922,1a + 9l2,2a) (3.13) 



+ yf i.9\p,a + 9\cc,p) + (5-^^)^ 9l2,29l2,cc 



2 
1 

+ [-2g^^g^^gi2,x + '7^^^7^' (2(?A;.,e - 9^.e,M^ ^2 ' 

-^3 = /b '^) (1 + '^^) v'^d^v, Ki = J^ip (x, t;) (1 + v'^) v'^d^v. 

Proof. See appendix B. □ 

The relations F^ = F^ = on G^ is to be arranged under a suitable choice 
of the distribution function on G^ = G^ x S. It is at this level that the 
Rendall method need to be modified. Assume that the distribution function 
if is such that 

T2a = i^aOnG'r^. (3.14) 

Then the reduced system i?2a = ^20 is equivalent to the following homoge- 
nous system of ODE on G^^ with unknown (F^, F^) 

9w'^% + {9^'^9i2,29w + 1^3/3,2) = 0, ,g 

g^pT% + (^l2^l2,2^74/3 + 1^4/3,2) r'' = ^ - ' 

Assumption f l3.14p is an integral system on in the sense that it is written 
explicitly as follows 

A^U + B^V = i,^, (3.16) 

where 

~ 8V2m2 ' ~ 8v^m2 (3.17) 

U = J^(p (x, f ) (1 + f ^) , y = (x, v) (1 + t;^) 

The determinant of the system ( 13.16^ is equal to *-~^g^4^ on G^^. So the 
solutions are given by the relations below 
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In brief has to be chosen in such a way that the following integral system 
(IS) holds for every x G 



(x, t;) (1 + v') = ^ {^P.Bs - ^,B,) , 

^ (x, t;) (1 + v^) v^d^v = ^ (^3^4 - ^4^3) . 

Assume F'' = on F. Then, in view of (I3.15p . F^ = on G^^. 



(IS) 



3.3. Relation F^ = on G^^. We seek for a combination of g^'^Rap, F^ 
and r\ that will provide an homogenous ODE on with unknown F^. 

3.3.1. Combination of g'^^Rap, F^ andT\. 
Proposition 8. On the following combination holds 

g'^^R^p - 2Vl - 2g^^g^2,2T^ = -^g^^ {N^p + M^p) , (3.21) 

where 

= —g^'^g^^g2\l {9l3fj.,a + 9fia,l3 — 9al3,fj.) „„s 
+9^^ {92PM + 92a,ip) - [^9^'^9l2,a + 9^^ (fl'M.a + fl-M^.A " 9aX,f,)] , ^ ' 

= [2g^^gu,x + g^^ (fi'/.e.A + gex,^^ - S'M.e)] [fi'''^ i9^tfs,a + fl-Ma,/? - fl-a/?,/.)] 

- i9^'^f (fl'12,/3 + 5'2/3,l) i9l2,a + 92a,l) 

- [9^^ {912,13 - 92p,l)] [^^^ (^12,a - ^2a,l)] 

- [9^^ i9ex,f3 + gei3,x - 9xpfi)\ [g^^ (s'^m," + 98a,,, - 9a^,&)\ ■ 

(3.23) 

Proof. It follows from relations (7.58 — 61) of by using (13.21) and (13.31) . 

□ 

3.3.2. Expression of g"'^ Tap. 
Proposition 9. On it holds that 

g'^^Tap = ^J^^ix,v) {v^fd^v 

+ ^^^^ (l - {^h,,)'^) V (x, v) v^v'd^v 

+ [{h^if (fi/iss)^ {(nhss)"^ - 2) + huh33] Is V (a;, v) {v^f d^v. 

(3.24) 

Proof. See appendix C. □ 
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In addition to assumption f l3.14p . assume 

g"^T^p = ^^/"^ + M,^) on G^^. (3.25) 

Then, in view of (]3.2ip . the reduced system Raj3 = Tap provides the following 
homogenous ODE on with unknown 

- 2Vl - 2(7! V2,2r' = 0. (3.26) 

The assumption 03.251) is a supplementary integral equation which is written 
explicitly as follows, Vx G 



^ (x, v) {v^ f d^v + telj^^ (1 - {nh,,Y) ^ (x, v) vWv 



{-9i2f 



= Ig'^P {Nafs + Map) . 

(IE) 

In view of fl3.26p . assuming = on F gives F^ = on G\,^. 

What we have just proved for the resolution of the constraints problem 
associated to the EV system can be summed up in the following main theo- 
rem. 

Theorem 1. Under the suitable integral assumptions {IS) and {IE) on the 
distribution function, there exists initial data for the reduced EV system such 
that the constraints F^ = are satisfied on U for the corresponding 
solution of the evolution problem associated to the EV system. 

Remark 6. It would be interesting to see whether the constraints integral 
equations can be avoided. One way of doing this is to work in temporal 
gauge and null moving frame (see But in the harmonic gauge case 

the issue may not be evident. Nevertheless we think that one could use an 
orthonormal frame in order to avoid that the metric appears (in an involved 
way) in the energy-momentum tensor (see [16j ). 



Appendix A: Proof of Proposition 6 

The first equality of (13. 4 p is provided in relation (7.42) of [10]. We handle 
the second one by using the expression (I2.28P of the energy-momentum tensor 
given in Proposition 5 to have 

T^J {x) = ^ j^^{x,v)v,, \g\'^ {-g'')~'^ IgH d\, (A.l) 

where Vij = f^. For i = j = 2 {A.l) reads 




(A.2) 
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From (13. ip we gain 



Thus 



P2 



V22 



-guq^ (1 + w^) on G\ 



1 



(A.3) 
(A.4) 

(A.5) 



1 2 



is now expressed on in terms of f via ( I2.20p to give 

The exact expression of (A^) is needed. Sphtting the quadratic form fl2.17p 
yields 

= (?22)^X2+ (^22)~^?2a^' 

+ (?33 - (922)"^ (?23)^) ' + (^22?34 " ?23?24) (?22 (?22?33 - (?23)^)) ' X'^ 

+ ((?22?33 - (?23)^) (?22?44 " (?24)^) - (?22?34 " ?23?24)^) ' (?22 (?22?33 - {923Y)) 

(A.6) 

By computing 'gij via tensorial transformation formulae we gain 

/ i^ll + 1^12 + i^22 1^11 - 1^22 I {913 + 923) I {9IA + 92a) \ 



(9. 



\ 



2i)^'-i I 4^tiz.i 4^911 4^922 2 

j9u - j922 j9ii - \9i2 + i5'22 2 (^'13 - fl'23) I (fi-u - 6-24) 

1 (fl'13 - fl'23) ^33 ^34 

2 (fl'14 — fl'24) ^34 5'44 



\ {913 + 923) 
2 (fl'14 + fl24) 



From ([23]) and ([S2D we get 

{9ij) 



( \9i2 








\ 





-|fll2 














fl33 


fl34 


I 





fl'34 


fl'44 / 



on G^. 



(A.7) 



Hence, in view of f l2.17p . we gain 

/ {-\9.2Y 



\ 










{^h^3) 








on G\ 



(A.8) 



33 



Thus, the inverse matrix of (A;^) reads 

/ (-|fll2)"^ 








{^3) 





\ 

-hi 



on G\ 



(A.9) 
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Simple calculation gives 

f^ = 2{g,^)-^ onG\ (A.IO) 

It is worth noting from (A.IO) that the condition < is equivalent to 
gi2 < on G^. We also have 

g^^ = OonG\ (A. 11) 
(A.5), (A.9), (A.IO) and (A.ll) imply 

= on G\ (A.12) 
It follows from (AA) and {A.12) that 

V22 = \{9i2f{l + v'Y onG\ (A.13) 

We now handle the term \g\^ i^^J^^) ^ \9\ ^ on G^. From (13. ip and (13. 2 p we 
gain 

det{gij) = -inguf on G\ (A.14) 

Thus 

\g\2 = Q\g^^\ = -Qg^^ on G^. (A. 15) 

Using (13. ip . (13. 2 p and tensorial transformation formulae we gain 

?22 = -^9l2, 92a = 0, gap = ^Kp OU 

Hence 

g = det (?ab) = -\9i2^'' on (A. 16) 

From (A.IO), (A. 15) and (A. 16) we gain 

W {-g'T^ Ian = I {912? on G\ (A.17) 

The insertion of (A.13) and (A.17) into (A. 2) gives the desired expression of 

T22 {x). 

Appendix B: Proof of Proposition 7 

(13. lip follows from relation (7.50) of by using gix = on G^. {x) 
is computed in the same manner as T22 {x). Actually, for i = 2 and j = a, 
(A.l) reads 

T2a{x) = -^ [ ^{x,v)v2a\9\H-9'T' lar^' dh. (B.l) 

2m^ Jb 

From (13. ip and (13. 2 p we gain 

Pa = q'9apw^ on G\ (B.2) 

(A.3) and (5.2) yields 

V2a = ^9l2 (1 + V^) {gaW + 9aiW^) OU G\ (B.3) 
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From ( I2I2OI) and {A.9) we gain 

= (-^12)5 {Qhs3yh!^ - huv^ , = ^{-gi2)^Q-'^ {h33)K*. 

(B.4) 

(5.3) and (5.4) imply 

t^2a = (1 + t^^) (-^12)^ ^ (ph33)~~^ - h34V^^ + Ki^-'^ (/^as)^ on G\ 

(B.5) 

Insertion of (A. 17) and (-B.5) into (5.1) gives the expression fl3.12p of {x). 

Appendix C: Proof of Proposition 9 
Tap {x) is handled in the same way as T22 and T2a- From {A.l) we have 

Tap{x) = ^ j^ip{x,v)vapW {-t^y'^W'^d^. (C.l) 

(5.2) implies 

fQ/3 = QaxgpixW^w^ on G\ (C.2) 
Insertion of (^4.17) and {C.2) into (C.l) gives 

Tap [x) = / (p (x, V) w^w^d'^v. 



B 



Thus 

g^^Tap (x) = ^^i^^ ^ (x, t;) {w' f + ^ (x, t;) u;^^;^^^^ 

(C.3) 

The insertion of (5.4) into (C.3) gives the expression (I3.24p of g'^^Tap (x) on 

Appendix D: Discussion on the integral constraints equations 
The integral system (IS). The integral system {IS) is written as follows 
{^xJ) = hi{x), {(p^,g) = h2{x) . (D.l) 
where (, ) denotes the scalar product in (5) and 

/(^) = (l+^;2)^3^ ^(^) = (l+^2)^4^ 

^1 = pSt^ (^3^4 - ^4^3) , h2 = ^ (^,^3 - ^3^4) . ^^-'^ 

As the distribution function must be non-negative, let us seek ip of the form 

(^) = ^ (a;, v) = [a (x) f {v) + h (x) g {v) + c {x)f , (D.3) 

where a (x) , h (x) and c (x) are unknown functions defined on C^. Expanding 
(5.3), we have 

ip = a^f + 52^2 ^ 2abfg + 2acf + 2bcg + (D.4) 
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where variables have been dropped for simphcity. Let ip be hke in {DA). It 
holds that 

/) = a\p, f) + h\g\ f) + 2ah{fg, f) + 2ac(/, /) + 2hc{g, f) + cP{l, /), 
{v., g) = a\p, g) + h^g\ g) + 2ah{fg, g) + 2ac(/, g) + 26c(^, g) + 0^(1, ^/). 

(D.5) 

Spherical coordinates will be used to calculate each of the following quantities 
that are needed. 

{p. f) = Ib [(1 + d'^^ id', f) = Ib [(1 + (1 + ^'d'^^ 

if 9, f) = L [(1 + v') v']' (1 + v') v^d\ if, f) = [(1 + d\ 
{9, f) = L (1 + ^'f ^'^'d\ (1, /) = L (1 + v') v'dh, 
{l,g) = (1 + v') v^d\ {g\ g) = [(1 + vf d\ 

{g,g)-^JAi^ + v')vfd'v. 

(D.6) 

As B is the open unit ball in R^, we set 

= r cos 6' cos A, f"^ = r cos6'sinA, f'' = rsin^, (D.7) 

with 

< r < 1, <e < A < 27r. (D.8) 

Define a domain P in as follows 

P = |(r, A) e RV < r < 1, - | < ^ < |, < A < 27r} . (D.9) 

Using the change of variables (-D.7), we gain 

(/^/) = /p/Vcos^dP, = Jpg'fr'cosedP, 

if 9, f) = /p cos ^dP, (/, /) = /p / cos edP, 

{gj) = Jpgfr'cosedP, (1,/) = Jpfr^ cos edP, (D.IO) 

(1'^) = J pgr^ COS 9dP, {g^,g) ^ fp g^r'^ cos 9dP, 

{g, g) = jp g'^r^ cos edP, 

where 

dP = drdOdX. 

After expansion and reduction we integrate the above quantities over P to 
obtain 



(D.ll) 



(/^ /) = {g\ f) = (/^ ^) = (5, /) = (1, /) = (1, 5) = {g\ g) = o, 

Prom {D.b) and (-D.11) it holds that 

647ra (a:)c(a;) 647rb {x)c{x) 

{^.J) = , {^.,g) = — ■ (D.12) 

Thus, if solves (D.l) if and only if 

647Ta{x)c{x) _ 647Tb {x)cix) _ 

105 105 ^^-^^^ 
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In sum if is given by 

ip^{v) =v{x,v) = [a{x) {l + v^)v^ + b{x) {1 + v^) + c {x)f , ... 
xeG\ v = {v\v\v^)eB, ^ > 

where 

a{x)c[x) = — ^4^, h[x)c[x) = — (D-15) 

The integral equation (/-E). The integral equation {IE) is written as 
follows 

^(^,, e) + /(^,, + 5(99,, s) = C, (D.16) 

where 

E{x)^^^, eiv) = {v'y 



(D.17) 



S (X) = (fi/l33)^ ((^^/i33)^ - 2) + /l44/i33" 

s{v) = {v')\ C{x) = \g'^^{N^^ + M^^). 
Using the expression of (p given in {DA) we gain 

{ip„ e) = a2(f , e) + e) + 2ah{fg, e) + 2ac(/, e) + 26c(5, e) + c\l, e). 

(D.18) 

As in the preceding paragraph, the above quantities are found to be 

(/2,e) = ||, (^^e) = |L, (l,e) = S, (/^, e) = (/, e) = e) = 0. 

(D.19) 

(D.18) and (D.19) imply 

, , Stt 9 Stt ,9 47r X 

Similarly it holds that 

{ip,, i) = a\f\ %)+h\g\ i)+2ab{fg, t)+2ac{f, t)+2bc{g, i)+c'{l, i). (D.21) 
Straightforward calculations as above give 

(/^ i) = (/, i) = (/, i) = (5, = (1, = 0, (/5, i) = ^ (D.22) 
(£1.21) and (D.22) give 

167ra6 

K,^) = ^■ (D.23) 
We are then left with calculating 

((^„ s) = a\f\ s) + 6^(52, s) + 2ah{fg, s) + 2ac(/, s) + 26c(^, s) + c\l, s). 

(D.24) 

By proceeding as above we get 

{f9,s) = {f,s)^{g,s)^Q, {f\s)^^, {9\s)^^^, {l,s) = ^^. 

(D.25) 
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(L).24) and (^.25) yield 

/ , Stt 9 87r,9 47r 9 x 

Prom {D.16) , {D.20) , {D.23) and (D.26), we see that the integral equation 
(IE) is equivalent to 

10 {3E + 3)0^ + 10 {3S + E)b'^ + 20Iab + 63{E + S) = ^C. (D.27) 

47r 

In view of (D.15), multiplying {D.27) by and rearranging , we gain 

110 250 (3^ (x) + S {hi {x)f + 110 250 (35 (x) + E [x)) (/i2 {x)f 
+220 500/ii {x) h2 (x) / {x) + 258 04871^ {E (x) + 5 {x}) [c (a;)]^ 
= 96 7 6 807rC (x) [c (a;)]^ 

(D.28) 

(L'.28) is an algebraic equation that can be solved under suitable assumptions 
to find c(,t). Doing so we deduce a{x) and b{x) thanks to (-D.15). Finally 
the distribution function if is obtained on and has the form {D.14). 

Acknowledgement. I am thankful to Professor Marcel Dossa who sug- 
gested this topic and guided me throughout the work. 
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